APS/123-QED 



O 

(N 
Oh! 

in 



Quasi-one-dimensional spin dynamics in LiV2 04: one-to-three dimensional crossover 

as a possible origin of heavy fermion state 
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Spin fluctuation in LiV2 04 is revisited by examining the earlier result of muon spin rota- 
tion/relaxation measurements. Instead of a relationship for the localized electron limit, one for 
itinerant electron systems between muon depolarization rate and spin fluctuation rate (vd) is em- 
ployed to re-analyze data, which reveals that vd varies linearly with temperature {ud oc T) over 
a range 10^-10^^ s~^ for 0.02 < T < 10^ K. Such a linear-T behavior as well as the magnitude 
of vd is fully consistent with that of the magnetic relaxation rate previously observed by inelastic 
neutron scattering (INS), demonstrating that /xSR and INS have a common time window over the 
fluctuation spectrum. The linear-T dependence of vd is understood as a specific feature predicted 
by a Hubbard model for intersecting one-dimensional (ID) chains. This quasi-lD character, which 
is co-existent with enhanced uniform susceptibility at low temperatures, supports the scenario of 
lD-to-3D crossover for the microscopic origin of heavy-fermion behavior in LiV204. 
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Heavy fermion (HF) behavior observed in a cubic vana- 
dium spinel, LiV204, has been in the spotlight of broad 
interest fii^ since it comprises one of remarkable examples 
in which only d-orbital electrons are relevant to the phe- 
nomenon. The formation of heavy quasiparticle (QP) 
state below a characteristic temperature, T* ~ 20 K, 
is suggested by large Sommerfeld coefficient (7 ~ 420 
mJ/molK^) and other bulk properties that are strikingly 
similar to typical /-electron HF compounds. Moreover, 
it has been inferred from the recent photoemission spec- 
troscopy that a peak of the density of states (DOS) just 
above the Fermi energy (Ep) develops for T < T*, which 
may correspond to the QP peak typically found for the 
/-electron systems 4 

While these observations seem to favor the Kondo 
mechanism established for the /-electron compounds as 
a common microscopic origin for the HF behavior in 
LiV204, approaches to support this scenario have been 
elusive. Earlier theoretical attempt to project the d- 
electron states (1.5 per V^'^"*" ion) onto the Kondo model 
by splitting them into two sub-bands by electronic cor- 
relationhad to introduce unusually large Kondo coupling 
{Jk ^ 10'^ K) to overcome competing effect of the Hund 
coupling. ^'^ Our /iSR study on single-crystalline sample 
provided evidence against the formation of a spin-singlet 
state, as it showed presence of staggered vanadium mo- 
ments at low temperatures far below T* that was inter- 
preted as the Kondo temperature {Tk) in this scenario,®'^ 
In the meantime, the importance of highly symmetric 
crystal structure and potential influence of geometrical 



frustration has been stressed by various authors, leading 
to a large variety of theoretical models^"— 

In the previous study using muon spin rotation and 
relaxation (/iSR), we have shown on a powder specimen 
of LiV204 that the observed ^SR signal consists of two 
components characterized by different response of depo- 
larization rate (A) to external magnetic field (/^o)'-- In 
particular, the signal with A showing least dependence 
on Hq (Ad, with a fractional yield / ~ 0.4) is mostly in- 
dependent of temperature below 10^ K, from which we 
suggested that the corresponding fluctuation rate derived 
from a general relation between A and 
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is also independent of temperature {i^d > 10^ s^^). 
In contrast, A associated with another signal (As, with 
1 — / ~ 0.6) is readily suppressed by Hq, which has been 
ascribed to slowly fluctuating local magnetic moments 
{vs ^ 10^-10^ Although the occurrence of such a 

phase separation has been confirmed by subsequent /iSR 
study on high-quality single-crystalline samples, the in- 
creased yield / (~ 0.8) strongly suggests that clarifying 
the origin of vd is essential to the understanding of elec- 
tronic state in LiV2 04fI 

Despite the fact that LiV2 04 is metallic, the use of 
Eq. ([T]) , which is valid in the limit of localized spins with 
h' determined by the local exchange interaction (J hv), 
is presumed to be justified by the presence of staggered 
vanadium moments suggested by broad /xSR linewidth 
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at low temperatures (which is also in line with pre- 
sumption of the Kondo scenario for T > T*)^ Besides 
this, Eq. (dl) may be regarded as a good approximation 
at high temperatures where the electronic state is sub- 
ject to strong damping by phonos. However, our recent 
/iSR study on another d-electron heavy fermion metal, 
Yi_3;Sca;Mn2, has demonstrated that more precaution 
should be taken for the interpretation of depolarization 
at lower temperatures 

Here, we revisit the spin dynamics of LiV204 on al- 
ternative basis that the staggered magnetic moments are 
carried by itinerant d electrons. The corresponding muon 
spin depolarization is described by a modified version of 
Eq. HD, 



(2) 



where the factors additional to Eq. ([T]) stems from elec- 
tronic density of state at the Fermi level with x be- 
ing the magnetic susceptibility, Na the Avogadro num- 
ber, and /is the Bohr magnetonJ^i2£l As a consequence, 
our reanalysis indicates that is linearly dependent 
on temperature {v oc T) with frequency ranging from 
10^-10^^ s~^, which accidentally serves to cancel the T- 
dependence of A. The distilled behavior of vjj turns 
out to be in excellent agreement with the previous re- 
sult of inelastic neutron scattering (INS)^ ^ strongly sug- 
gesting that both ^SR and INS have been observing a 
common phenomenon. Moreover, the linear-T depen- 
dence of spin fluctuation, which is commonly observed in 
Yi_2;Sca;Mn2, is understood as a property characteristic 
to the spin-spin correlation of the intersecting Hubbard 
chains that simulate the pyrochlore lattice.— This implies 
the crucial role of t2g orbitals as one-dimensional (ID) 
chains that are under a strong geometrical constraint of 
pyrochlore lattice structure, and further suggests the di- 
mensional crossover due to coupling between these chains 
as the primary origin of the heavy-fermion statCi^ 

Since the previous study on single- crystalline sam- 
ple indicates that muon depolarization is dominated by 
a component showing fast fluctuation {vd), we focus 
on the behavior of vd- The evaluation of using 
Eq. ([2]) requires additional information on magnetic sus- 
ceptibility (x) shown in Fig. [TJ which has been ob- 
tained on the powder sample used for the previous ^SR 
measurement.^ Considering the behavior of % observed 
in single-crystalline sample^^i^^ we attribute the diver- 
gent behavior with T to unknown paramagnetic 
impurities and decompose data into two parts. 
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where 6 is the Weiss temperature. Curve fits using 
Eq. © yields a = 0.387(5) emu-K/mol, 6 = 74(2) K, and 
b = 0.0131(2) emu-K/mol, in which the behavior of xv 
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FIG. 1: (Color online) Magnetic susceptibility (x) as a func- 
tion of temperature in the powder specimen used for /iSR 
measurement in Refii. Solid curve shows the best fit using a 
form described in the text. 



is in good agreement with that of single crystals. Hence- 
forth we use X — Xv for the input of Eq. ^ . Assuming 
that the Curie term comes from free V'^'^+ spins (l.S/is), 
the fractional yield of impurity phase estimated by 6/T 
is 1.6% of the total volume. This is far smaller than that 
of primary /iSR signals (either / or 1 — / with / ~ 0.4), 
indicating that the paramagnetism of the impurity phase 
is irrelevant to the interpretation of /iSR data. 

Another important quantity in Eq. ([2]) is the hyper- 
fine parameter, 5^. The muon Knight shift measure- 
ments in both powder and single-crystalline samples yield 
(5^(£)) ~ 0.5 ± 0.2 GRz/ fiB that corresponds to the sig- 
nal with vd^^ Apart from the large error due to broad 
linewidth, it is in good agreement with calculated value, 
^niD) = 0.143 GHz//iB, for muons that occupy a site 
in the center of cyclic vanadium hexamer (as inferred 
from the observed /iSR linewidth due to nuclear mag- 
netic moments at high temperatures) and are subject to 
the magnetic dipolar fields from vanadium ionsj^ Since 
the experimental value of hyperfine parameter has large 
uncertainty due to broad linewidth, we use the calculated 
value for the evaluation of i^d below. 

Fig. [5] shows muon depolarization rate (Xd) under zero 
external field, deduced by curve fits using a sum of two 
exponential damping^ 

P,(<) = /exp(-ADi) + (1 - /) cxp(-A5<). 

Although the data are scarce particularly at higher tem- 
peratures. Ad is only weakly dependent on temperature 
with a tendency of gradual increase with decreasing tem- 
perature, which is qualitatively similar to the behavior 
of Xv (see the inset of Fig. These features are re- 
markably similar to that observed in Yi_a;ScxMn2 with 
greater Sc content {x > 0.07).^^* 

The evaluation of from A^i is straightforward by 
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FIG. 2: (Color online) Muon depolarization rate in LiV204 
reproduced using data published in Ref|^ for the component 
showing fast fluctuation {fo)- Inset: reproduction of Fig. [1] 
with logarithmic scale for temperature. 
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FIG. 3: (Color online) Spin fluctuation rate (y) as a function 
of temperature in LiV204. Thin solid line shows linear T 
dependence {v ocT). Inelastic neutron scattering (INS) data 
are quoted for comparison, where the dashed curve shows the 
linewidth (Fg-Fo)//! at Qc = 0.64 A"^ with Fo ~ 1.5 meV.— 
(See text for more detail.) 



using Eq. ([2]) in a modified form, 
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where Hq = for the present condition of zero external 
field. The result is plotted in Fig. [3] together with data of 
INS,^^ where one can observe that I'd falls on a straight 
line representing linear relation to temperature (vd T) 
over a T range of three decades. This is again strikingly 
similar to the behavior of spin fluctuation rate observed 
in Yi_^Sc:^:Mn2 with x > 0.07.^^ 

According to a numerical simulation of quarter-filled 
Hubbard chains (or rings) by path integral quantum 
Monte Carlo method, the spin dynamics is predicted to 
exhibit a linearly T-dependent relaxation rate in the spin- 
spin correlation. 



= 1^0 + c(q) • T, 



(5) 



where vq comes from residual ferromagnetic spin corre- 
lation (corresponding wave vector q = 0).^> It is also 
inferred from the simulation that Eq. ([5]) depends only 
weakly on q and other parameters such as U /t (where 
U and t are the on-site Coulomb energy and the trans- 
fer integral between nearest neighboring vanadium sites 
in the Hubbard model, respectively). The observed be- 
havior of I'D is perfectly in line with the above predic- 
tion with i/Q ~ 0. As shown in Fig. 3, the magnitude 
of vd as well as its linear-T dependence is in excellent 
agreement with the spin relaxation rate observed in INS 
at q ~ Qc = 0.64 A^^ (around which a broad peak is 
observed below T*) after subtracting the residual term, 
r,(T ^ 0) = To ~ 1.5 meV^ The difference between 
vo (mSR) and Tq (INS) may be attributed to that in the 
sensitive q range of observation, as the q-integrated Fg 



(for 0.6 < q < 1.3 A^^) shows tendency of smaller Fqj^ 
It must be remembered, however, that is subject to 
certain ambiguity due to the broadness of linewidth 
observed in the muon Knight shift measurement. 

While the comparison between i^d and F^ demon- 
strates that both /iSR and INS probe nearly identical 
part of the spin fluctuation spectrum in LiV204, it also 
confirms that ^Li-NMR has a completely different time 
window of observationi^l^— Considering estimated hy- 
perfine field on ^Li nuclei {Am ~ 0.258 T/^s in Refl25l, or 
0.171 T/fJ-B in Ref,f26l), the corresponding parameter 6^ 
in Eq. (g]) is evaluated to be 7LiAhf ^ 17.8-26.9 MUz/fiB 
(where 'ju = 27r x 16.546 MHz/T). Since the spin- lattice 
relaxation rate (1/Ti = A) is reported to vary over a 
range of 10^-10^ s~^ and to follow the Korringa relation 
(A/T ~ 2.0-2.25 s-^K-i) for T < T* ~ 20 K^i^-^ the 
corresponding spin fluctuation rate evaluated by Eq. ^ 
levels off to take a constant value, I'u ^ 10^'^ s^^ 
(T ^ T*). Interestingly, the estimation of vu using re- 
lationship between the electronic specific heat coefficient 
and N{Ef), 7 = (271^ / 3) klN{EF)N a = 0.420 J/mol 
yields 
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hN{Ei. 
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= 1.7 X 10" s~^ 
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This is in excellent agreement with the experimental ob- 
servation, and strongly suggests that l/Ti in ^Li-NMR 
is dominated by the Pauli paramagnetism. It should be 
stressed that the above frequency is beyond the limit of 
observation for ^SR under the present condition, because 
the corresponding A would be less than 10~^ MHz. 

A theoretical model by Fujimoto presumes the quasi- 
ID character of the t2g bands associated with the py- 
rochlore lattice (consisting of intersecting chains of t2g 



4 



orbitals) as an essential basis for the description of elec- 
tronic state in LiV204, because the hybridization be- 
tween the ID bands will be strongly suppressed owing 
to the geometrical configuration.^'^ It is interesting to 
note that such quasi- ID bands may have singularities in 
DOS, as in an analogous case of Al5-type intermetal- 
lic compounds (e.g., VaSi) where it is predicted that 
N{E) oc 1/{E - £'c)^/^ with being the position of 
the d bandit The model incorporates the hybridization 
as a perturbation to the ID Hubbard bands, which yields 
an energy scale (T*) that characterizes the dimensional 
crossover from ID to 3D as the Fermi-liquid state devel- 
ops with decreasing temperature below T* . The calcu- 
lated 7 taking account of the latter as the leading correc- 
tion to the self-energy yields a large value, consistent with 
the experimentally observed ones. The progression of hy- 
bridization also induces the enhancement of the 3D-like 
spin correlation that would appear as the enhancement of 
uniform susceptibility, while the spin fluctuation is domi- 
nated by the staggered component of ID Hubbard chains. 
The increase in x (= Xv) with decreasing temperature 
shown in Fig. [1] is understood as the manifestation of 
such a dimensional crossover, as it perfectly follows the 
prediction by the relevant theory The quasi- ID char- 
acter of the low-energy spin fluctuation preserved below 
T* coexisting with the enhanced x, which is commonly 
observed in Yi_a;Scj;Mn2, comprises strong evidence of 
such a scenario. 

It would be honest to mention that the muon depo- 
larization rate (A, from which v is deduced) does not 
provide any criterion in itself to decide the credibility of 
the present interpretation based on itinerant model com- 
pared with that based on the local spin model adopted 
in our previous report. However, the itinerant model has 
substantial merit as it provides a coherent understanding 
of spin fluctuation in metallic LiV204 without resorting 



to the splitting of t2g band for a part of vanadium ions 
to bear "local" spins, where such a band splitting is yet 
to be demonstrated experimentally. 

Finally, we briefly mention the signal component with 
slower spin fluctuation observed by /xSRi^ii This com- 
ponent is characterized by relatively small hyperfine pa- 
rameter ((5^(s) = 23-34 MHz//iB), from which fluctuation 
rate I's is evaluated by Eq. (|H) to yield 10^-10'' s~^ at low 
temperatures (T < 1 K). Considering the fact that the 
fractional yield of this component decreases drastically in 
the single-crystalline samples (1 — / ~ 0.2), one may be 
inclined to associate this component to some unknown 
impurity phase. However, there is no other experimental 
evidence for such a second phase having different chem- 
ical form with the relative yield as high as ^20%. Thus 
more careful investigation would be needed to clarify the 
origin of this slow fluctuation. 

In conclusion, we have shown that reanalysis of muon 
depolarization rate in LiV204 using the model of spin de- 
polarization for the itinerant electron systems yields spin 
fluctuation rate that is linearly dependent on tempera- 
ture. This result turns out to be fully consistent with that 
of the previous inelastic neutron scattering experiment, 
providing a coherent understanding of the spin dynam- 
ics by the theoretical model of intersecting ID-Hubbard 
chains that simulates pyrochlore lattice. The persistent 
quasi- ID spin dynamics coexists with the enhanced uni- 
form susceptibility at lower temperature (T < T*), which 
is common to another d-electron heavy-fermion system 
Yi_2;Sca;Mn2. These observations strongly suggest that 
geometrically constrained t2g band is the primary stage 
for the formation of heavy quasiparticles with the ID- 
to-3D dimensional crossover as a possible mechanism of 
effective mass enhancement. 
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